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We have made a topological study of cosmic microwave background (CMB) polarization maps by simulating the 
AMiBA experiment results. A ACDM CMB sky is adopted to make mock interferometric observations designed 
for the AMiBA experiment. CMB polarization fields are reconstructed from the AMiBA mock visibility data using 
the maximum entropy method. We have also considered effects of Galactic foregrounds on the CMB polarization 
fields. The genus statistic is calculated from the simulated Q and U polarization maps, where Q and U are Stokes 
parameters. Our study shows that the Galactic foreground emission, even at low Galactic latitude, is expected to 
have small effects on the CMB polarization field. Increasing survey area and integration time is essential to detect 
non-Gaussian signals of cosmological origin through genus measurement. 
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I. INTRODUCTION 

The cosmic microwave background (CMB) polariza- 
tion can provide useful information in determining cos- 
mological parameters that are only weakly constrained 
by the CMB temperature anisotropy alone, such as 
the epoch of reionization or the presence of tensor 
perturbations. The Array for Microwave Background 
Anisotropy (AMiBA; Lo ct al. 2001), an interferomet- 
ric array of 19 elements, will have full polarization capa- 
bilities in order to probe the polarization properties of 
the CMB. In this paper, we simulate the AMiBA polar- 
ization experiment and study the topological properties 
of the AMiBA mock maps through the genus statistic 
of Q and U Stokes parameter fields. We also investigate 
the effect of polarized Galactic foreground emission on 
the CMB polarization fields. 

II. Simulating the AMiBA Polarization Ex- 
periment 

(a) Observational Strategy 

An efficient observational strategy is essential to im- 
age the CMB polarization field. Although Subrah- 
manyan (2001) suggests drift scanning that is more de- 
sirable in reducing the environmental response and the 
cross-talk, we adopt simple observational strategy that 
the co-mounted array can be rotated and shifted to the 
adjacent points after some data acquisition. We assume 
that the primary beam A(x) of an elemental aperture 
with diameter 0.3 m is Gaussian, normalized to unity 
at the peak and with the FWHM of 44', the bandwidth 
of the dual polarization receivers is 20 GHz centered at 
90 GHz, and the system temperature is 75 K. When 
0.3 m apertures are used, the AMiBA will be sensitive 
to CMB polarization over the range 700 < £ < 2000 
with signal to noise ratio of about 4 at ^ — 700 and 2 



at ^ - 1150 in 24 hours (Lo ct al. 2001). 

We assume that 19 elements are hexagonally close- 
packed on a platform, as shown in Figure la, where 
the distances between adjacent aperture centers are 35 
cm, giving 171 baselines with the minimum spacing of 
35 cm. This configuration gives a good sampling of 
the uw-plane, and is useful for reconstructing the CMB 
polarization maps from the visibility data. 

In a single pointing observation, the visibilities are 
assumed to be measured for 2 hours at a set of u speci- 
fied by the close-packed configuration, then the instru- 
ment is rotated by 5° about an axis through the center 
of the aperture plane to obtain a different set of u with 
the same baseline lengths but with different orienta- 
tions. In this way the uw-plane will be well-covered, 
allowing for good imaging. The number of rotations 
are chosen to be 12 times per pointing for better sam- 
pling of the uw-plane. The uv coverage after a 24 hour 
observation is shown in Figure 16. Figures Ic and Id 
show the window function in the I space and the beam 
pattern arising from our simulation, respectively. For 
large I ^ 60, we can relate spherical harmonic multi- 
pole index i with a baseline coordinate u by £ = 27rM 
with a good approximation. 

The resolution in the u space is limited by the area 
of the sky that we have surveyed. This is equal to the 
size of the primary beam in a single pointing observa- 
tion. We can increase the resolution of the map as well 
as the survey area by mosaicking several contiguous 
pointing observations (see White et al. 1999 for more 
details). Mosaicking does not increase the range of u. 
It simply enhances our resolution by allowing us to fol- 
low more periods of a given wave, which is analogous 
to the Fraunhofer diffraction through many holes. In 
our simulated observation, we perform a mosaic survey 
of 9 X 9 pointings with spacing FWHM/2. 
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Fig. 1. — (a) A hexagonally close-packed configuration for AMiBA experiment, 
(6) the uv coverage a 24 hour observation, (c) the window function in the £ space, 

and (d) the beam pattern for a given set of uv coverage in a 2? 5 x 2? 5 field. A 
thick circle in the middle shows the FWHM (44') of the primary beam A(x). 
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(b) Simulated Observations 

The visibility function sampled at a pointing po- 
sition Xp on the sky is the Fourier transform of the 
true sky brightness /(x) weighted by the primary beam 
A{x - Xp), i.e., 



y(u, 



,)7(x)e^ 



(1) 



(see, e.g., Ng 2001; Hobson, Lasenby, & Jones 1995). 
The size of the primary beam A(x) determines the area 
of the sky that is viewed and hence the size of the map, 
while the maximum spacing determines the resolution. 

In order to simulate an observation of CMB polar- 
ization fields, we follow the prescription given in Seljak 
(1997) and Zaldarriaga & Scljak (1997). Wc use the 
CMBFAST package (Scljak & Zaldarriaga 1996; Zal- 
darriaga, Seljak, & Bertschinger 1998) to get the polar- 
ization power spectrum (Cf) and the cross-correlation 
between temperature and polarization (Cf) in a flat 
ACDM model with Ha = 0.6, = 0.4, h = 0.6, and 
flb/i^ = 0.0125 (Ratra et al. 1997) in which only the 
E-mode (scalar mode) polarization exists and the B- 
mode vanishes. To generate small patches of Stokes 
parameter Q and U fields (5° x 5°) on the sky we use 
small-scale limit approximation by Zaldarriaga & Sel- 
jak (1997) 
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Q{x)= j dPu[E{u) cos20u-S(u) sin2</)u]e" 

U{x) = j dPu[E{u) sin 2</)u + B{u) cos 2(/)u]e""^" '', 

(2) 

where £'(u) is the Fourier component in u space of E- 
polarization field, 4>u is the direction angle of the two- 
dimensional vector u, and B{u) = 0. Figure 2a shows a 
realization of a pure CMB Q field in the ACDM model 
in a 5° X 5° patch on the sky. 

According to equation (1), these fields are then mul- 
tiplied by the primary beam and Fourier transformed 
to give a regular array of visibilities. An AMiBA ob- 
servation is simulated by sampling the regular array 
at the required points in the uii-plane specified by the 
AMiBA observational strategy (see Fig. 16). 

The AMiBA instrument noise on the data is simu- 
lated by adding a random complex number to each vis- 
ibility whose the real and imaginary parts arc drawn 
from a Gaussian distribution with the variance of the 
noise predicted in a real observation. Here we use 
the sensitivity per baseline per polarization defined as 
(Wrobel & Walker 1999; Ng 2001) 
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"Hsria Aphys y/2Al'Tacc ' 

iwidth, Tacc is the correla 
lation time, Aphys is the physical area of an elemental 



where Av is bandwidth, Tacc is the correlator accumu- 



aperturc, and rjs (rja) is system (aperture) efficiency. 
Figure 2b shows a dirty image directly FFTed from the 
single pointing mock visibility data. Since the dirty im- 
age contains biased field information, it cannot be used 
for topology analyses. Thus we need to reconstruct 
the original and unbiased CMB field from the observed 
noisy visibility data. 

(c) Image Reconstruction Using the MgLximum 
Entropy Method 

Prom the mock visibility data, CMB polarization 

fields can be reconstructed using the maximum entropy 
method (MEM), in which an extended emission field 
can be more effectively restored compared to the com- 
mon CLEAN algorithm (see Narayan & Nityananda 
1986 for a review, and see Conwell & Evans 1985; 
Chae & Yun 1994 for applications). The MEM is also 
more appropriate for reconstructing a non-Gaussian 
field, compared to the Wiener-filtering method (see, 
e.g., Maisinger et al. 1998). We use the methods of 
Cornwell & Evans (1985) for general the MEM algo- 
rithm, and of Cornwell (1988) for MEM mosaicking. 
However, the standard MEM method contains a loga- 
rithmic term that is inapplicable to images that have 
both positive and negative values, such as fluctuations 
in the CMB. Therefore, we consider the image to be the 
difference between two positive additive distributions, 
namely hi = fi — gi, and the entropy becomes 



S{h,m) = J2 



ijji — 2mj — hi In 



2m,- 



(4) 



where hi is the predicted map, nii is the default model 
map (Maisinger, Hobson, & Lasenby 1997; Jones et al. 
1998), and = [/if + 4mf]i/2. 

The attraction of MEM is that it deflnes the best 
image as the solution of a simple optimization problem 
obtained by maximizing the entropy while fitting to the 
data. Thus extra information can simply be added as 
constraints in the optimization. In our application of 
mosaicking, is defined as 
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Here V denotes the predicted visibility, and ^ is the 

variance in the visibility of the rth visibility sample at 
the pth pointing (see Cornwell 1988 for more details). 
Maximizing 

J = S - ax^, (6) 

where a is a Lagrange multiplier, gives a solution for 
h that maximizes the entropy and fits to the data. We 
find the most appropriate image iteratively using the 
Newton- Raphson method, i.e.. 



Ah = (-VV J)-^ • VJ. 



(7) 
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Q Polarization 



(a) Rire CMB (d) Galaxy 1 f 1 ,b)= f 90,20 ) deg,- ja gm. 




(b) Dirty: Sii^le Poind t^. Noise (e) GaJaiy Q; theta_c = 0-5 deg 




Fig. 2. — Examples of mock CMB Q polarization fields, Galactic foregrounds, 
and MEM reconstructed fields. All numbers given in colorbars have an unit of 
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Here we choose the defauh model map as m = 5 fiK, 
constant everywhere, and the Lagrange multipher as 
a = 100. These values are insensitive to the final so- 
lution. The mosaicked images reconstructed by MEM 
are very close to the original CMB polarization maps 
(compare Fig. 2c with Fig. 2a). Figure 3 compares 
the CMB Q and U fields along a row in the original 
map (dotted curves) with those in the MEM map re- 
constructed from 9x9 mosaic data (solid curves). 

(d) Galactic Foregrounds 

Wc investigate the effects of polarized Galactic fore- 
ground emission on the CMB polarization fields. Two 
possible sources of the polarized Galactic emission are 
synchrotron and spinning dust emissions, which are 
known to induce at least 10 % polarization fraction at 
frequencies below 90 GHz (Lubin & Smoot 1981). Vi- 
brational dust and free-free emissions are expected to 
be unpolarized. We use a simple toy model of Kogut & 
Hinshaw (2000) to make foreground polarization maps. 
We assume here that mean fractional polarization is 10 
%, the polarized Galactic emission is proportional to 
the unpolarized intensity, and the effect of synchrotron 
emission is negligible at 90 GHz. We can model the 
Stokes Q and U components as 

g = /cos(27)7, C/ = /sin(27)7, (8) 

where f{0) is the fractional polarization, assumed to 
vary across the sky with (/) =0.1, and j{d) is the po- 
larization angle which is randomly realized with coher- 
ence angle 6c = 0?5 (see §2 of Kogut & Hinshaw 2000 
for more details). We use a two-component dust model 
of Finkbeiner, Davis, & Schlegel (1998) to predict 90 
GHz vibrational dust intensity, and a model of Kogut 
et al. (1996) for adding the free-free emission. The 100 
/Ltm dust map of Schlegel, Finkbeiner, & Davis (1998) 
is used as a basic template map. Figures 2d, 2e, and 
2/ show the 90 GHz Galactic intensity at 6 = 20°, its 



Q polarized contribution, and the MEM reconstructed 
CMB Q polarized field contaminated by the Galactic 
polarized foreground, respectively. 

III. Topology of the Simulated Maps 

The CMB polarization field, like the temperature 
anisotropy, can be used for an observational test of the 
Gaussianity of the primordial density fluctuation, and 
thus provide an important constraint on the inflation- 
ary models. We use the two-dimensional genus statistic 
introduced by Gott ct al. (1990) as a quantitative mea- 
sure for the topology of the CMB polarization fields. 
The genus is defined as the number of hot spots minus 
the number of cold spots, or equivalently, 

where k is the signed curvature of the iso-temperature 

contours C. Since the genus curve as a function of 
the temperature threshold level has a form of g{iy) 

oc i'e~^ for a Gaussian random-phase field, non- 

Gaussianity of a field can be detected from deviations 
of the genus curve from this relation. For example, 
non-Gaussianity can cause shift or asymmetry of the 
genus curve. 

We present genus curves as a function of the area- 
fraction threshold level i^a, which is defined to be the 
temperature threshold level at which the correspond- 
ing iso-temperature contours enclose a fraction of the 
survey area equal to that of a Gaussian field (Gott 
et al. 1990). All measured values are averages over 
15 mock AMiBA surveys. We use the 4 FWHM x 4 
FWHM region in each mock survey field, where the 
MEM reconstruction gives more confident results. Fig- 
ures 4a and 46 shows averaged histograms and genus 
for foreground- free (solid curve). Galactic foreground- 
added (open circles), and Galactic foreground- added 
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Fig. 4. — Histograms and genus measured from CMB polarization Q maps. 



and MEM reconstructed (9x9 mosaic observation, 
filled circles) CMB Q fields. The last case is a result of 
a 24 hour integration per pointing located at Galactic 
latitude h = 20°. Note that the solid curve in Figure 
46 shows the functional form expected for a random- 
phase Gaussian field, Ave~^ and has been fitted to 
the foreground-free genus result by adjusting A. The 
result from a 24 hour integration of a region centered 
at 6 = 35° is shown in Figure 4c. When the integration 
time is only one hour per pointing, we obtain genus 
curves shown in Figure Ad. 

Table 1 lists skewness {{Q/a)^), genus-related statis- 
tics A (amplitude), Aiy (shift parameter), and Ag 
(asymmetry parameter) measured from CMB Q fields; 
see Peebles (1993) for skewness and Park et al. (2001) 
for definition and measurement of genus-related statis- 
tics, respectively. 

The results show that the polarized Galactic emis- 
sions, even at low latitude, do not have significant effect 
on the CMB polarization fields. Although the polar- 
ized Galactic foreground emissions at 6 = 20° cause 
slight asymmetry, and reduce genus amplitude (Fig. 46 
and Table 1), these effects are smaller than the statisti- 
cal fluctuations due to the sample variance of the pure 
CMB map. The open and filled circles in each panel 
in Figure 4 show that the original polarization fields 
can be well restored from the AMiBA mosaic observa- 
tional data by MEM image reconstruction technique if 



the integration time is 24 hours per pointing. However, 
if the integration time is one hour, instrumental noise 
increases the amplitude of the genus curve, and the 
presence of non-Gaussian signals due to the Galactic 
emission is buried. 

IV. Conclusions 

We have studied the topology of CMB polariza- 
tion maps by simulating the AMiBA experiment and 
by measuring the skewness and genus statistics. The 
MEM turns out to be very useful in reconstructing the 
CMB polarization Q and U fields from the interfcro- 
mctric visibility data. Our study shows that although 
the polarized Galactic foreground emissions at low lat- 
itude (6 = 20°) decrease the amplitude of the genus 
curve and cause asymmetry of the curve, these effects 
arc not significant compared with the sample variance 
in our 2? 9 x 2? 9 surveys. At higher Galactic latitudes, 
the effects of the Galactic emissions i s shown to be neg- 
ligible, and the intrinsic non-Gaussianity of the CMB 
polarization is easier to detect. 

The major source of statistical variance of the genus 
statistic is the sample variance induced by our small 
survey area (4 FWHM x 4 FWHM = 8.6 deg.^). This 
can be reduced not by observing many small different 
patches of the sky, but by increasing the area of the 
survey at each location. Therefore, increasing the sur- 
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Table 1. Skewness and Genus Statistics Measured from Mock CMB Q Fields 





Q Pol. 




A 


A^/ 


A.g 




CMB 


-0.08 ±0.11 


52.5 ±4.4 


-0.006 ±0.035 


0.04 ±0.10 


6 = 20° 


+ Galaxy 


-0.06 ±0.16 


47.9 ±4.7 


-0.005 ±0.047 


0.04 ±0.12 




+ Galaxy (observed)''' 


-0.05 ±0.15 


48.6 ±4.6 


-0.014 ±0.044 


0.03 ±0.12 


6 = 35° 


+ Galaxy 


-0.07±0.13 


51.6 ±3.9 


-0.009 ±0.040 


0.05 ±0.11 




+Galaxy (observed)^ 


-0.06 ±0.12 


52.9 ±4.2 


-0.008 ±0.037 


0.02 ±0.09 



9x9 mosaic observation with 24 hours integration time per pointing. 



vey area and integration time at a given pointing is 
essential to accurately estimate the effects of Galac- 
tic foregrounds on the polarization maps and to detect 
non-Gaussian signals of cosmological origin. 

Wc acknowledge valuable discussions with Kin- Wang 
Ng, Uros Seljak and Ravi Subrahmanyan. This work 
is supported by the BK21 program of the Korean Gov- 
ernment. 
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